Wigner symmetry, Large N c and Renormalized One Boson Exchange Potentials. 
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Wigner symmetry in Nuclear Physics provides a unique example of a non-perturbative long 
distance symmetry, a symmetry strongly broken at short distances. We analyse the consequences of 
such a concept within the framework of One Boson Exchange Potentials in NN scattering and keeping 
the leading JV C contributions. Phenomenologically successful relations between singlet So and triplet 
3 Si scattering phase shifts are provided in the entire elastic region. We establish symmetry breaking 
relations among non-central phase shifts which are successfully fulfilled by even-L partial waves and 
strongly violated by odd-L partial waves, in full agreement with large N c requirements. 
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I. INTRODUCTION 

Symmetries have traditionally been very useful in nu- 
clear physics partly because the force at the hadronic 
level is not well known at short distances [H, @, In 
some cases, like isospin, chiral or heavy quark symmetry, 
the invariance can directly be traced from the fundamen- 
tal QCD Lagrangean and formulated in terms of the un- 
derlying quark and gluonic degrees of freedom. In some 
other cases the connection is less straightforward. Many 
years ago Wigner and Hund proposed [3, Q extending 
the spin and isospin SUs(2) g) SUj(2) symmetry into the 
larger ST/ (4) group where the nucleon-spin states p f, 
p J,, n |, n I correspond to the fundamental representa- 
tion, and hence providing a supermultiplet structure of 
nuclear energy levels as well as new selection rules for 
nuclear transitions and response functions Q. The cor- 
responding SU (4) mass formula was found to be at least 
as good as the well known Weizsacker one 0,1- Spin- 
orbit interaction of the shell model obviously violate the 
symmetry, and indeed a breakdown of SU(A) has been 
reported for heavier nuclei 0] while nuclear matter has 
been addressed in [To| . Double binding energy differences 
have been shown to be a useful test of the symmetry [ll[ . 
Recently, inequalities for light nuclei based on SU (4) and 
Euclidean path integrals have been derived by neglecting 
all but S-wave interactions fl2l ]. 

Despite its relative success along the years, ST/ (4) sym- 
metry has been treated as an accidental one within the 
traditional approach to Nuclear Physics and guessing its 
origin from QCD has been a subject of some interest 
in the last decade. Indeed, attempts to justify SC/(4) 
spin-flavour symmetry from a more fundamental level 
have been carried out along several lines. Based on the 
limit of large number of colors N c of QCD [lU, it 
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was shown (lH [ill ] that if the nucleon momentum scales 
as p ~ N®, the nuclear potentials scale either as N c or 
l/N c , depending upon the particular spin-isospin chan- 
nel, which shows that the NN force could be determined 
with 1/iV 2 relative accuracy. It was found that the lead- 
ing potential would be SU (4) symmetric if the tensor 
force was neglected in addition, a plausible assumption 
for light nuclei where S-waves dominate. Although these 
estimates arc conducted directly in terms of quarks and 
gluons, quark- hadron duality allows to reformulate these 
results in terms of purely hadronic degrees of freedom, 
providing a rationale for the One-Boson-Exchange (OBE) 
potential models [13], and the internal consistency of 
Two- PJ] and Multiple Boson Exchanges [H, The 
analysis of sizes of volume integrals of phenomenologi- 
cally successful potentials confirms the large N c expec- 
tations [2l| . The large size of scattering lengths was re- 
garded as a fingerprint of the SU (4) symmetry within 
an Effective Field Theory (EFT) viewpoint [22| using 
the Power Divergent Subtracted (PDS) scheme; singlet 
and triplet renormalized couplings coincide at the natu- 
ral renormalization scale fj, ~ m n 3> l/cv s , 1/ott, with a s 
and at the scattering lengths, and a contact interaction 
makes sense in such a scaling regime. Resonance satura- 
tion based on the elimination of exchanged mesons in the 
OBE Bonn potential [l_7| at very low energies was also 
shown to reproduce the EFT approach and to agree nu- 
merically with the Wigner symmetry expectations [23[. 
According to Refs. [H [2| HI S3, HI QCD might be close 
to a point where the effective theory had an SU (4) sym- 
metry at zero energy as well as discrete scale invariance if 
the pion mass was larger than its physical value, around 
TOtt ~ 200MeV. This nice idea might be confirmed by re- 
cent fully dynamical lattice QCD determinations of the 
scattering lengths [29[ and quenched lattice QCD evalu- 
ations of NN potentials [3(J EH where indeed unphysical 
pion masses are probed. 

While the proclaimed symmetry holds in a range where 
scale invariance sets in and EFT methods based on con- 
tact interactions can be applied [13, [H| , is not obvious 
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what are the implications for the lightest NN system it- 
self for finite energies and for physical pion masses. In 
particular, the scale dependence of the contact interac- 
tion is modified when the finite range of the long distance 
potential is taken into account. To be specific, low en- 
ergy NN scattering is dominated by S-waves in different 
(S, T) channels where spin and isospin are interchanged, 
(1,0)= 1 5 ^ (0, 1)= 3 Si. Wigner SU{4) symmetry pre- 
dicts identical interactions in both 1 Sq and 3 Si channels. 
The above mentioned identity of the 1 Sq and 3 Si po- 
tentials holds also in the large iV c -expansion [l5|, [l6| , so 
we take advantage of this fact by using the leading N c - 
OBE potentials which simplifies the discussion to a large 
extent as we discuss in Sect. [Til In contrast, the corre- 
sponding phase shifts from Partial Wave Analyses [32j 
are very different at all energies. We are thus confronted 
with an intriguing puzzle since it is not obvious at all in 
what sense should the symmetry be interpreted for the 
NN system; it would be difficult to understand otherwise 
the successes of SU (4) for light nuclei. A second puzzle 
arises from an embarrasing cohabitation of conflicts and 
agreements between large- N c studies and Wigner sym- 
metry. Despite the initial claim [l5j a more complete 
analysis [l6j could only justify the Wigner symmetry in 
even-L partial waves while for odd-L a violation of the 
symmetry was expected. However, doing so required ne- 
glecting the tensor force, which according to the Wigner 
symmetry should vanish, but it is a leading contribu- 
tion to the potential in the large N c limit. Thus, while 
some pieces of the NN potential (such as e.g. spin-orbit) 
are suppressed in both schemes, some others are not si- 
multaneously small. These conflicts between the time- 
honoured SU (4) Wigner symmetry and the QCD based 
large- N c expansion for odd-L channels require an expla- 
nation and naturally pose the question on the validity of 
either framework. 

In the present work we analyze both puzzles by intro- 
ducing the concept of a long distance symmetry firstly 
to understand the meaning of Wigner symmetry in those 
cases where its validity complies with large N c expecta- 
tions. This is a case where we expect the symmetry to be 
more robust. Once this is done, it is pertinent to diluci- 
date the validity of the symmetry in those cases where a 
possible conflict with the large N c expansion arises. Our 
discussion is tightly linked to the coordinate s pac e renor- 
malization discussed in previous works [HI, |34| . This 
approach while borrowing the physical motivation from 
EFT theories provides a quantum mechanical framework 
where the emphasis is placed on the non-perturbative as- 
pects of the NN problem, a playground where the stan- 
dard EFT viewpoint has encountered notorious difficul- 
ties. The method is reviewed in Sect. IIIII for complete- 
ness. We find that for S-waves the Wigner symmetry 
holds in a much wider range than the applicability of a 
contact interaction suggests if the finite range of the in- 
teraction is incorpored. As a byproduct we provide in 
Sect. ITVl quantitative predictions; the seemingly indepen- 
dent triplet and singlet S-waves phase shifts correspond- 



ing to iso- vector and iso-scalar states respectively for the 
np system are shown to be neatly intertwined in the en- 
tire elastic region. A similar correlation can also be estab- 
lished between the 1 S'o virtual state and the 3 Si deuteron 
bound state. Actually, we show how the symmetry may 
be visualized for large scattering lengths due to the on- 
set of scale invariance. Symmetry breaking due to in- 
clusion of further counter-terms, tensor interaction and 
spin-orbit interaction are discussed in Sect. [Vj We show 
how a sum rule for supermultiplet phase shifts splitting 
due to spin-orbit and tensor interactions is well fulfilled 
for non-central L-even waves, and strongly violated in L- 
odd waves where a Serber-like symmetry holds instead. 
This pattern of S[/(4)-symmetry breaking complies to 
the large N c expectations, a somewhat unexpected con- 
clusion. Finally, in Sect. fVTl we provide our main conclu- 
sions and outlook for further work. 



II. OBE POTENTIALS, LARGE N c AND 
WIGNER SYMMETRY 

Our starting point is the field theoretical OBE model 
of the NN interaction [ItJ which includes all mesons with 
masses below the nucleon mass, i.e., ir, <r(600),7y, p(770) 
and w(782). For the purpose of discussing SU(A) Wigner 
symmetry within the OBE framework (see Appendix 
for a brief overview) we will deal here with S-waves only, 
neglecting for the moment the S-D wave mixing stem- 
ming from the tensor force as required by Wigner sym- 
metry. Our results of Sect. IIVI and estimates in Sect. IV Bl 
will provide the a posteriori justification of this simplify- 
ing assumption. Non-central waves and the role of spin- 
orbit as well as tensor force will be treated in Sect. IV CI 
as S/7(4)-breaking perturbations. For the S-waves the 
NN potential reads 



V = Vc + tW c + o-Vs + to-Ws , 



(1) 



where r = t\ • ti = 2T(T + 1) — 3 and a = o~\ ■ o-i — 
2S(S + 1) — 3 and Pauli principle requires (— ) S + T + L = 
— 1. Thus, for the spin singlet and spin triplet 3 Si 
states we get 



V s = V c + W c - 3V S - 3W S , 
V t = Vc-3Wc + Vs-SWs, 



(2) 
(3) 



To simplify the discussion we will discard terms in the 
potential which are phenomenologically small. Actually, 
according to Refs. [la, [l6| in the leading large N c one 
has Vc ~ W s ~ N c while V s ~ Wc ~ l/N c . In terms of 
meson exchanges (see also Ref. [l8|) one has the contri- 
butions 



V.(r)=V t (r) = - 



9l N N m le m - r 
l6irM% r 



4-7T r 
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pNN 
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0{N, 



(4) 
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where g a NN is a scalar type coupling, g^NN a pseudo- 
scalar derivative coupling, g^NN is a vector coupling and 
fpNN the tensor derivative coupling (see [l?} for nota- 
tion). Here, the scheme proposed in [35| of neglecting 
both energy and nonlocal corrections is realized explic- 
itly. In principle the large iV c limit contains infinitely 
many multi-meson exchanges which decay exponentially 
with the sum of the exchanged meson masses. However, 
NN scattering in the elastic region below pion production 
threshold probes CM momenta p < p max — 400 MeV. 
Given the fact that l/m w = 0.25fm <C {/p m ax = 0.5fm 
we expect heavier meson scales to be irrelevant, an in 
particular to and p themselves, are expected to be at 
most marginally important . Note that, in any case, 
when to w = m p the redundant combination g^AW — 
f 2 NN m 2 /(2M^) appears, indicating a further source of 
cancellation between p and to in this channel. Moreover, 
since the leading contributions to the potential are ~ N c 
and the subleading ones are ~ 1/N C) the neglected terms 
are of relative 1/iVj order, so we might expect an a priori 
~ 10%, accuracy. 

The coincidence between 1 Sq and 3 «Si potentials com- 
plies to the Wigner SU (4) symmetry which we review 
for completeness in Appendix [X] for the two-nucleon sys- 
tem. Modern high quality potentials [IB] describing ac- 
curately NN scattering below pion production threshold 
show some traces of the symmetry for distances above 
1.4 — 1.8fm. Quenched lattice QCD evaluations of NN 
potentials for m^/nip ~ 0.6 [30l.l3lT| yield also similar 1 Sq 
and 3 Si potentials for r > 1.4fm. Thus, at first sight one 
may conclude that Wigner symmetry holds when OPE 
dominates, and thus has a limited range of applicability. 
An important result of the present investigation, which 
will be elaborated along the paper, is that this is not 
necessarily so, provided the relevant scales of symmetry 
breaking are properly isolated with the help of renormal- 
ization ideas. 

Let us analyze the consequences of the symmetry, 
Eq. ([3]), within the standard approach to OBE poten- 
tials. The scattering phase-shift So(p) is computed by 
solving the (S-wave) Schrodingcr equation in r-space 



U ;(r)+M N V(r) u p (r) 
u p (r) 



p 2 u p (r), (5) 
sin (pr + S (p)) 



sin S (p) 



, (6) 



with a regular boundary condition at the origin u p (0) = 
0. Moreover, for a short range potential such as the one 
in Eq. ^ one also has the Effective Range Expansion 
(ERE) 



pcot S {p) 



1 



■ r P 



(7) 



where the scattering length, ao, is defined by the asymp- 
totic behavior of the zero energy wave function as 



Uo(r) 



r 

1 

a 



and the effective range, ro, is given by 

2 

I I r \ 

ro = 2 



dr 



I-—) -u (r) 2 
a 



(8) 



(9) 



In the usual approach [17|,[37| everything is obtained from 
the potential assumed to be valid for < r < oo. We note 
incidentally that the Wigner symmetry relation, Eq. Q , 
holds at all distances 2 . In addition, due to the unnatu- 
rally large NN So scattering length (oe s ~ — 23fm), any 
change in the potential V — ► V + AV has a dramatic 
effect on cto, since one obtains 



Aao = cyqMn 



AV{r)u (r) 2 dr, 



(10) 



and thus the potential parameters must be fine tuned, 
and in particular the short distance physics. As it was 
discussed in Refs. [38l l39j this short distance sensitivity 
is unnatural as long as the OBE potential does not truely 
represent a fundamental NN force at short distances. In- 
deed, the sensitivity manifests itself as tight constraints 
for the potential parameters when the 1 Sq phase shift 
is fitted resulting in incompatible values of the coupling 
constants as obtained from other sources as NN scatter- 
ing. Of course, there is nothing wrong in the need of a 
fine tuning as this is a unavoidable consequence of the 
large scattering length; the relevant point is whether this 
should be driven by a potential which will not be realistic 
at short distances. 

In any case, in the traditional approach to NN po- 
tentials we are confronted with a paradox; on the one 
hand the symmetry seems to suggest that somewhere the 
phase shifts should coincide, while on the other hand a 
fine tuning is required because of the large scattering 
lengths. In the standard approach, if V s (r) = V*(r) then 
5 s (k) = St(k) and thus a s = at, as one naturally expects. 
A straightforward explanation, of course, is to admit that 
the symmetry is strongly violated. This would make dif- 
ficult to understand how can SU (4) work at all for light 
nuclei if the simpler two nucleoli system does not show 
manifestly the symmetry. 

Before presenting our solution to this dilemma in the 
next section, let us note that a good condition for an ap- 
proximate symmetry is that it be stable under symmetry 
breaking, otherwise a tiny perturbation V s {r) — Vt(r) = 
AV(r) would yield a large change, and this is pre- 
cisely the bizarre situation we are bound to evolve be- 
cause of the large scattering lengths. This suggests a clue 



This of course does not exclude explicit and leading N c uncorre- 
cted multiple pion exchanges at, i.e. background non-resonant 
contributions in tttt or wp scattering. We expect them not to be 
dominant once a, p and 10 are included. 



2 In practice, strong form factors are included mimicking the finite 
nucleon size and reducing the short distance repulsion of the 
potential, but the regular boundary condition is always kept. 
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to the problem, namely we should provide a framework 
where the highly potential-sensitive scattering length be- 
comes a variable independent of the potential. More 
generally, we want to avoid the logical conclusion that 
a symmetry of the potential is a symmetry of the S- 
matrix 3 . As we will explain below the puzzle may be 
overcome by the concept of long distance symmetry; a 
symmetry which is only broken at short distances by a 
suitable boundary condition. 



III. UNIVERSALITY RELATIONS AND 
RENORMALIZATION 

We cut the Gordian knot by appealing to renormaliza- 
tion in coordinate space [33l |34| . As we will show this 
enables to disentangle short and long distances in a way 
that the symmetry is kept at all non- vanishing distances. 
The main idea is to fix the scattering length indepen- 
dently of the potential by means of a suitable short dis- 
tance boundary condition. As a result the undesirable 
dependence of observables on the potential is reduced 
at short distances, precisely the region where a determi- 
nation of the NN force in terms of hadronic degrees of 
freedom becomes less reliable. 

Let us review in the case of S-waves the renor- 
malization procedure in coordinate space pursued else- 
where [33L l34j and which will prove particularly suitable 
in the sequel. This is fully equivalent to introduce one 
counter-term in the cut-off Lippmann-Schwinger equa- 
tion in momentum space (see Ref. [4l[ for a detailed dis- 
cussion on the connection). The superposition principle 
of boundary conditions implies, 



u k (r) = Uk, c (r) + k cot 5 u k , s (r) , 



(11) 



with Uk,c{f) — ► cos(kr) and Uk, s (f) ~ * sin(kr)/k for r 
00. At zero energy, k — > 0, and 8o(k) — > — agk yields 



uo(r) = u .c(r) uo, s (r) 

a 



(12) 



with uo,c(f) — ► 1 and uo. s (r) — » r for r — ► 00. Combining 
the zero and finite energy wave functions we get 



[u' k {r)u (r) - u' (r)u k (r)] 



u k (r)u {r)dr . 



(13) 



where r c is a short distance cut-off radius which will be 
removed at the end. To calculate the contribution from 



3 This situation resembles the case of anomalies in Quantum Field 
Theory where the parallel statement would be that a symmetry 
of the Lagrangean becomes a symmetry of the S- matrix, a conclu- 
sion which may be invalidated by the impossibility of preserving 
the symmetry by the necessary regularization of loop integrals. 
The present case is a bit more subtle as it corresponds to the 
case of finite but ambiguous theories (see e.g. Ref. |40|) 



the term at infinity we use the long distance behavior, 
Eq. ©. The integral and the boundary term at infinity 
yield two canceling delta functions. This corresponds to 
take 



u k (r)u p (r)dr = 



7T<5(fc — p) 

2 sin 2 S (k) 



(14) 



as can be readily seen. We are thus left with the bound- 
ary term at short distances, taking the limit r c — > we 
get 



lim [u' k (r c )u (r c ) - u' (r c )u k (r c )} = . 



(15) 



Note that the regular solution u k (r c ) = Uo(r c ) = is a 
particular choice for r c — 0. Writing out the orthogonal- 
ity condition via the superposition principle at finite and 
zero energies, Eq. (ITT1) and Eq. (TT2"]) respectively, one gets 







dr 



uo,c(r) uo, s (r) 

a 



x u k , c (r) + k cot 5 (k) u k , s (r) 
Expanding the integrand and defining 

/>oo 

A(k) = / dr uo, c (r)uk, c (r) , 



(16) 



B(k) 
C{k) 
V(k) 



dr u ,s(r)u k , c (r) , 
dr u ,c(r)u kiC (r) , 



dr u , s (r)w M (r) , 



we get the explicit formula 

a a A{k) +B{k) 



kcotSo(k) 



a C(k)+V(k) 



(17) 



(18) 



The functions A, B, C and T> are even functions of k 
which depend only on the potential. Note that the de- 
pendence of the phase-shift on the scattering length is 
wholly explicit] cot <5o is a bilinear rational mapping of 
ato- Further, using Eq. (fl"2")) , one gets the effective range 



A 



B 



C 



a a\ ' 



where 



.4 
B 
C 



2 / dr(l - u 2 J , 
Jo 

poo 

-4 / dr(r - u , c uo, s ) , 
Jo 



= 2 



dr(r 2 - u 2 ) . 



(19) 

(20) 
(21) 
(22) 



depend on the potential parameters only. Again, the 
interesting thing is that all explicit dependence on the 
scattering length ao is displayed by Eq. (Ti~9l ) . 
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We turn now to discuss the case of a bound state cor- 
responding to the case of negative energy E = —^ 2 /M 
where 7 is the wave number. The wave function behaves 
asymptotically as 

u 7 (r)^A s e^ r 7 (23) 
and is chosen to fulfill the normalization condition 



u 7 (r) 2 dr = 1 . 



(24) 



In principle, such a state would be unrelated to the scat- 
tering solutions. An explicit relation may be determined 
from the orthogonality condition, which applied in par- 
ticular to the zero energy state yields 







dr 



1 



uoA r ) u , a {r) 



0() 



uJr) . 



(25) 



This generates a correlation between the scattering 
length, olq and the bound state wave number, 7, 



ao(7) 



/ °° drM 7 (r)u 0;S (r) 
/ °° dru 7 (r)u , c (r) 



(26) 



We remind that the two independent zero energy solu- 
tions, u c (r) and uq s (r) depend only on the potential. 

A trivial realization of the conditions discussed above is 
given by the case where there is no potential, U(r) = 0. 
Hence, the general solution for a positive energy state 
E = k 2 /M is given by 



Uk(r) = cot So(k) sin(fcr) + cos(fcr) . 



(27) 



and using the low energy limit condition 8o(k) — > — a^k 
we obtain 



u o(r) = 1 

a Q 



(28) 



Orthogonality between zero and finite energy states 
yields after evaluating the integrals 



kcotSo{k) = — — . 

a 



(29) 



and as a consequence the effective range vanishes r = 0, 
in accordance to the fact that the range of the potential 
is zero. For a negative energy state E = —j 2 /M the 
normalized bound state is 



u 7 (r) = A s e~ 



A s = l/v^y. 



(30) 



Orthogonality between the zero energy and the bound 
state, again, yields the correlation 



a = I/7. 



(31) 



In the appendix iBl we illustrate further the procedure in 
the case of weak potentials for which a form of perturba- 
tion theory may be applied for the case of weak potentials 
but arbitrary scattering lengths. 



Wigner correlation rJl/Un) 



4 r— 

3.5 -S 

3 - 

2.5 - 

2 - 

1.5 - 

I - 

0.5 - 

— 
-0.4 



71-exch 
(jl + CT)-exch 
r(th) =2.695 
r (exp) = 2.770 
r t (th) = 1.628 
r t (exp)= 1.753 




-0.2 



0.2 



0.4 



0.6 



l/otfl [fin 



FIG. 1: The Wigner correlation for the effective range ro 
(in fm) of a np S-wave as a function of an arbitrary inverse 
scattering length ao in the case of the OPE and OPE+cr po- 
tentials. The parabolic shape is determined by a unique long 
distance potential. The points in the solid curve correspond 
to the two different values of the effective range r B in the sin- 
glet 1 So and r t in the triplet 3 Si channels when the scattering 
length is taken to be a s = — 23.74fm and a t = 5.42fm respec- 
tively. Experimental points are also shown for comparison. 



Before going further we should ponder on the need to 
take the limit r c — > 0, which corresponds to eliminating 
the cut-off. We note that the potential, V(r), is used 
at all distances both in the standard approach, which 
involves the regular solution only, and the renormaliza- 
tion approach, which requires the regular as well as the 
irregular solution. However, the sensitivity to the short 
distance behaviour of the potential is quite different; the 
standard approach displays much stronger dependence 
while the renormalization approach is fairly independent 
on the hardly accessible short distance region, a feature 
that becomes evident perturbatively (see e.g. Eq. (|B6|) h 
This is in fact the key property that allows to eliminate 
the cut-off in the renormalization approach. Thus, re- 
moving the cut-off does not mean that the OBE potential 
is believed to hold all the way down to the origin. 

The procedure carried out before is described in purely 
quantum mechanical terms, but it can be mapped onto 
field theoretical terminology; it is equivalent to the 
method of introducing one counter-term in the cut-off 
Lippmann-Schwinger equation in momentum space (4ll . 
l42j . Moreover, Eq. (fT2"|) represents the corresponding 
renormalization condition, which is chosen to be on- 
shell at zero energy. In the case of the bound state 
the corresponding renormalization condition is given by 
Eq. (|23[) at negative energy. Imposing more than one 
renormalization condition, i.e. introducing more than 
one counter-term and removing the cut-off presents some 
subtleties which have been discussed in Refs. [H, l4ll |. 
We will analyze below this issue in the present context 
(see Sect. ESty - 
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IV. CENTRAL PHASES AND THE DEUTERON 



Potential Parameters 



To proceed further we fix the potential parameters, 
keeping in mind that the leading N c nature of the poten- 
tial embodies some systematic 1/N 2 uncertainties. Of 
course, while we will use relations which are compatible 
with large N c scaling, the numerical values can only be 
fixed phenomenologically. The main point is that besides 
the er-meson mass (see below), we may choose quite nat- 
ural values for the masses and couplings unlike the usual 
OBE potentials As was discussed at the end of 

Sect. |TT] the standard approach suffers from tight con- 
straints reflecting the unnatural short distance sensitiv- 
ity. In this regard, let us note that, as emphasized in 
Refs. [HI, [3^1, it is a virtue of the renormalization view- 
point which we are applying here to the OBE potential, 
that the unwanted short distance sensitivity is largely 
removed, allowing for a determination of the potential 
parameters using independent sources. For definiteness 
we take g^NN = 13.1 and g^NN = 10.1, quite close to the 
Goldberger-Treiman values for a and 7r, g a NN = Mn / 
and g-nNN = gAM^/f^ respectively. We also take the 
SU(3) value g^NN — 3gpNN — g<f>NN which on the ba- 
sis of the OZI rule, g<f>NN = 0, Sakurai's universality 
g P NN = 5 P 7r7r/2 and the KSFR relation 2g 2 f 2 = m 2 



Liiiversal function A 



I VCl'S: UlllUltill B 



yields g u 



NN 



N c m p /(2V2U) = 8.8. The rho tensor 



coupling is taken to be J p nn — V^M^g^j^N /m p = 15.5 
which cancels the vector meson contributions in the po- 
tential and yields k p = f P NN I gpNN = 5.5 a quite reason- 
able result [17| 4 - Note that 1/N C effects include not only 
other mesons but also finite width effects of a and p since 
for large N c one has stable mesons, T a ,T p ~ 1/N C . For 
the masses we take = 140MeV and m w = 783MeV. 
This fixes all parameters except m a (actually the real 
part) which we identify with the lightest J PC = ++ me- 
son / (600). According to the recent analysis based on 



Roy equations m a - iT a /2 = 441+£° - i272tj 2 MeV [43j. 
A fit to the pn data of Ref. (36[ in the 1 Sq channel yields 
m a = 510(l)MeV, where the error is statistical. The fit- 
ted mass value differs by about 10% from the location of 
the real part of the resonance, in harmony with the ex- 
pected 1 /N 2 corrections 5 . Although a more quantitative 
estimate of the large N c corrections to the potentials pa- 
rameters would be very useful, for the present purposes 



As shown in previous work 38, 39] the net vector meson exchange 
contribution corresponding to the combined repulsive coupling 
91nn ~ fpNN m2 p/ 2M N ( referred ther e simply as g* NN ) cannot 
be pinned down accurately from a fit to the 1 Sq phase shift 
being compatible with zero within errors. This is due to the short 
distance insensitivity embodied by the renormalization approach. 
Actually, our estimate of the cr-mass as a pole in the second 
Riemann sheet for tttt scattering for large Nc [3p| yields the value 
m a ~ 507MeV. 




50 100 150 200 250 S 
PjMeV] 

Universal function C 



50 100 150 200 250 300 350 400 450 
PcnJMeV] 

Universal function D 




FIG. 2: The universal functions A, B C and T> defined by 
Eqs. (|17[) in appropriate length units as a function of the CM 
momentum p (in MeV) . These functions depend on the poten- 
tial V„(r) = Vt(r) only but are independent of the scattering 
length. 



of discussing Wigner symmetry on the light of large N c it 
is more than sufficient. Thus, we make no attempt here 
to make any systematic expansion. 



B. Low energy parameters and phase shifts 



Clearly, in the traditional approach if we have V s (r) = 
Vt(r) and impose the regular boundary condition, 
M s (0) = Uf(0) = 0, the only possible solution is a s = a t , 
r s = rt and S s (p) = 6t(p). However, in the renormaliza- 
tion approach we allow different short distance boundary 
conditions <(0+)/w s (0+) ^ u' t (0 + )/u t (0 + ) 6 , and hence 
we may have a s ^ a t . Note that this corresponds to a 
breaking of the symmetry at short distances and hence 
postulating its validity at long distances. The previous 
equations imply straight away the following expressions 



6 The limit from above, it(0 + ) = lim r _^q+ u(r c ) is really necessary 
to pick both the irregular and irregular solutions. If one starts 
exactly from the origin the only possible solution is the regular 
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for the effective ranges in the singlet and triplet channels, 



A 



D 



0:. s 



c 



n = 



, B C 

°4 



a, 



(32) 



As already mentioned, the remarkable aspect of these 
two equations is the fact that the coefficients A,B,C are 
identical both in the triplet as well as in the singlet chan- 
nels as long as V s (r) — Vt(r), thus the only difference 
resides in the numerical values of the scattering lengths 
a s and a t - Numerically we get (everything in fm ) 



r = 1.3081 

ao 

= 1.5089fm (ao = 
= 0.6458fm (a = 

, = 2.4567- 5 ' 5284 



4.5477 5.1926 



a s ) (exp.2.770fm) 
at) (exp.l.753fm) 
5.7398 



a 

2.6989fm (a 
1.5221fm (a 



(tt + a) 



a s ) (exp.2.770fm) 
at) (exp.l.753fm) 



(33) 



where the corresponding numerical values when the ex- 
perimental a s — — 23.74fm and at — 5.42fm as well as 
the experimental values for the effective ranges have also 
been added. More generally, for any fixed potential the 
correlation of tq on 1/ao is a parabola which we plot 
in Fig. Q] for the OPE and OPE+o\ This dependence 
is universal to all S-waves having the same potential and 
from this viewpoint there is nothing in this curve making 
unnaturally large scattering lengths particularly different 
from smaller ones. The present analysis, however, does 
not shed any light on the origin of the large size of the a's 
nor how a s and at are interrelated 7 . In any case, as we 
see from Fig. [TJ the experimental values fall strikingly 
almost on top of the curve, pointing towards a correct 
interpretation of the underlying symmetry. 

We turn next to the phase shifts. According to Eq. (Tig)) 
they are given in terms of the universal functions A., B 
C and V defined by Eqs. (fT7|) and presented in Fig. [3] 
in appropriate length units as a function of the CM mo- 
mentum p in MeV for completeness. As we see, these 
functions are smooth. From them the corresponding sin- 



Nijmegen 
71 -exchange 
G-exchange 
(7l+G)-exchange 




Pcm [MeV] 




50 100 150 200 250 300 350 400 450 
Pcm [MeV] 

FIG. 3: Phase shifts (in degrees) for the fitted 1 So (top panel) 
and predicted 3 Si (bottom panel) channels as a function of 
the CM momentum (in MeV). In both cases the potential 
is the same, V s (r) = Vt(r), while the only difference is in 
the scattering lengths in the singlet channel a s = — 23.74fm 
and in the triplet channel a t = 5.42fm corresponding to a 
different short distance boundary condition. We also plot the 
cases with only lcr-exchange and l7r-exchange for comparison. 
Data from 



glet and triplet phase shifts are obtained by 



k cot S s 
k cot St 



a s A(k)+B(k) 



a s C(k) 
a t A(k) 



V(k) 
-B{k) 



a t C(k)+V(k) 



(34) 



7 This is in fact a price we pay for the built-in short distance 
insensitivity. We note, however, that after Refs. 24, 25, 26, 27, 
|2§H both scattering lengths might coincide for a pion mass around 
m w ~ 200MeV. As a consequence, QCD might be close to a 
point where the effective theory had a standard SU(A) symmetry 
at zero energy. Actually, in Ref. |25| the similarity between 1 So 
and 3 Si phase shifts can be seen. This scenario would turn the 
long distance symmetry we propose for the physical pion mass 
into a standard symmetry for such an unphysical value of the 
pion mass. 



respectively. When the experimental scattering lengths 
a s = — 23.74fm and at — 5.42fm are taken we can fit the 
singlet 1 5'o channel and predict the triplet 3 Si channel. 

The result is shown in Fig. [3] and as we see the agree- 
ment is remarkably good taking into account that we 
have neglected the tensor force and the a priori l/N% 
systematic corrections to the potential. Note that the 
identity of the singlet and triplet potentials is not suffi- 
cient; the simple OPE fulfills this property but does not 
explain the neither phase-shifts. Actually, it shows that 
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both failures are correlated 8 . 



C. Renormalization and scale invariance 

It is interesting to analyze our results on the light of 
Refs. [TUHH, where a square well potential, PDS and 
sharp momentum cut-off were used respectively to model 
the short distance contact interactions arising when all 
exchanged particles are integrated out. Here we are inter- 
ested in the dependence on the arbitrary renormalization 
scale separating the contact and the extended particle ex- 
change interaction since they are not independent of each 
other; by keeping this scale dependence we may enter the 
interaction region where, as we will show now, the sym- 
metry can be visualized. We appeal to the coordinate 
space version of the renormalization group [U, [44[ (for a 
momentum space version see Ref . [45[ ) , where the version 
of the Callan-Zymanzik equation for potential scattering 
reads 

Rc' {R) = c (i?)(l - c Q (R)) + MR 2 V(R) , (35) 

where Cq(R) — Ru' Q (R) /uq(R) is a suitable combination 
of the short distance boundary condition and we have 
chosen for simplicity to work at zero energy 9 . The above 
equation provides the evolution of the boundary condi- 
tion as a function of the distance R (the renormalization 
scale ) in order to have a fixed scattering amplitude (see 
Ref. 34[ for a thorough discussion). Clearly, at long dis- 
tances r 3> I/jtItt the potential becomes negligible and 
the equation is scale invariant, only broken by the renor- 
malization condition which fixes the value of cq at some 
scale 10 . In fact, the solution of the above equation is 
given in terms of the scattering length ao in the infrared, 
R — ► oo , co(-R) — > ao/R- On the other hand, if the 
scattering length is large we also have an intermediate 
regime with clear scale separation and 



co(R) 



R — ao 



-1 . 



\jm v < R < a , (36) 



indicating the onset of scale invariance [34j. This is in 
agreement with the PDS argument of Refs. [22j if the 
identification n ~ 1/R is done. Eventually the infrared 
stable fixed point cq — > 1 will be achieved. Note, however, 
that c s (R) ~ Ct(R) in a much wider range, particularly in 
the scaling violating region where the potential acts. In 
the more conventional language of wave functions the sit- 
uation corresponds to a case where both wave functions 




FIG. 4: Zero energy S-waveradial functions for the singlet 1 So 
and triplet 3 Si channels as a function of distance (in fm). The 
normalization is such that ti ig — * 1 — r/a s and u 3g 1 — > 
1 — r/at with a s = — 23.74fm and at = 5.42fm the singlet and 
triplet scattering lengths respectively. The potentials gener- 
ating these wave functions are the same V„(r) = Vt(r). 



u s (r) ~ ut(r) for r <C a Sl at- The situation is illustrated 
in Fig. 3] where the similarity in the range below lfm can 
clearly be seen, and does not differ much from the solu- 
tion uo tC (r) entering the superposition principle, Eq. (fT^]) 
and corresponding to the limit ao — > ±oo. Note that the 
symmetry can be visualized within the range of the po- 
tential only when the scattering length is large because 
there exists the scaling regime <C r <C «0i but the 

long distance correlations between the two S-wave chan- 
nels due to the identity of potentials hold regardless of 
the unnatural size of the scattering lengths. 



D. Virtual and bound states 

It is of course tempting to analyze the kind of features 
for the deuteron that may be obtained from this simpli- 
fied picture where the tensor force is neglected from the 
start. The deuteron is determined by integrating in the 
Schrodinger equation with negative energy E = —^\jM 
with 7 d = 0.2316fm _1 the wave number and imposing 
the long distance boundary condition, Eq. (|23|) . We also 
compute the matter radius 



2 1 
T = — 

m 4 



r 2 u d {rY 



(37) 



8 The reason why OPE fails at much lower energies in the 1 5o 
channel than in the 3 Si channel is due to a stronger short dis- 
tance sensitivity of the channel with larger scattering length. 

9 The orthogonality conditions discussed above correspond to take 
c p (R) -* c (R) for R -> 0. 

10 In appendix \U\ we analyze a case where the dilatation symmetry 
of a 1/r 2 potential must necessarily be broken by a renormaliza- 
tion condition. 



and the M.mi matrix element 



A S M M \ = / dru d (r)u i So (r) . 
Jo 



(38) 



which correspond to the dominant magnetic contribution 
to neutron capture process np — > jd in the range of ther- 
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mal neutrons (<■ 



KeV) in stars 11 . 

-l 



For the experimen- 
tal ld = 0.2316fm _i we get A s = 0.8643fm- 1/2 (exp. 
0.8846(9)fm~ 1/2 ) and r m = 1.9138fm (exp. 1.9754(9)fm 
and M M i = 4.0464fm) (exp. 3.979fm). As mentioned 
above, orthogonality between the bound state and the 
zero energy state yields an explicit correlation between 
the triplet scattering length, at and the deuteron wave 
number, 7, 



J °° dru~,(r)uo, s (r) 



/„ dru 7 (r)uo !C (r) 



-y=id 



(39) 



Since the two independent zero energy solutions, uo iC (r) 
and ito, s (r) depend only on the potential and hence are 
identical for the S-wave components of the singlet and 
triplet channels, this correlation is a consequence of the 
Wigner symmetry as well as long as we take V s (r) = 
Vt(r). Note that taken as a function of the scattering 
length, the expression 



M(y,a ) 



dr w 7 (r)uo(r) , 



(40) 



yields both the orthogonality relation as well as M. 

M(jd,a t ) = 0, 
M(-y d ,a s ) = M M i ■ 



Ml 



(41) 



Actually the dependence on the inverse scattering length 
is a straight line which we show in Fig. [5l As we see both 
conditions are very well fulfilled. Similarly to the previ- 
ous case, the orthogonality between finite energy states 
and the deuteron corresponds to the magnetic contribu- 
tion to the photodisintegration of the deuteron. The re- 
sult, however does not differ much from the potential- 
less theory, and so we will not discuss it any further. 
For the experimental 7^ = 0.2316fm we get a t — 
5.32fm. This value improves over the simple formula 
at = I/7 = 4.31fm obtained from the case without po- 
tential, or the single OPE case where at = 4.60fm. It is 
worth stressing that the same relation above yields the 
virtual state, a purely exponentially growing wave func- 
tion, u v {r) — » e +7vr , in the singlet channel, yielding for 
a s = -23.74fm -the value -y v = 0.042fm _1 . In other 
words, the function 00(7) fulfills «o(7d) = ct t and simul- 
taneously Q!o(7v) = ck s . Numerically we get 



a (-0.042fm" 
a (0.2265fnT 



= -23.74fm : 
= 5.42fm. 



(42) 
(43) 



11 In this normalization the total cross section is given by 

°m("P -» yd) = ira([i p - fi n ) 2 VB/2E(B/M N )~/M 2 M1 

where Eis the neutron energy, and /i p and fi n the proton and 
neutron magnetic moments in units of the nuclear magneton, 
MJV = e/(2M p ). We are neglecting meson exchange currents in 
the calculation of A4m±. 



Ml matrix element 



(ji + o)-exch 
7C-exch 
no potential 
Ml'V) = 4.047 
Ml exp (CL) = 3.979 
Ml(a,) = 




-0.1 



-0.05 



0.05 

l/oc [fm" 



0.1 



0.15 



0.2 



FIG. 5: The matrix element M(y, a) for 7 = 7^ as a func- 
tion of the inverse scattering length ao. This function should 
fulfill the properties A4(yd,& s ) ~ Ml mi (the neutron capture 
Ml matrix element) and M("fd, at) ~ (the orthogonality re- 
lation between the deuteron and the zero energy triplet state. 
The experimental values are highlighted as points. 



In the region below lfm the virtual state u v (r) and the 
deuteron bound state Ud(r) look very much alike the cor- 
responding singlet and triplet zero energy wave functions 
respectively (see Fig. 0]). Thus, u ^s (r) ~ u v (r) and 
u 3 Sl (r) ~ Ud(r) are consequences the closeness of the 
poles to the real axis, either in the second or first Rie- 
mann sheets respectively. However, u Q x s(r) ~ u 3 s (r) 
and u v (r) ~ u d (r) are further consequences of the iden- 
tity of the potentials V s (r) — V t (r) . 



V. SYMMETRY BREAKING 
A. Symmetry breaking with two counter-terms 

An essential ingredient of the present analysis is the 
requirement of orthogonality between different energy 
states, which ultimately reflects the self-adjoint charac- 
ter of the Hamiltonian. This implies that, for the Yukawa 
like potentials we arc dealing with, the only way to pa- 
rameterize the unknown information at short distances 
is by allowing, besides the regular solution, the irregular 
one and fixing the appropriate combination by impos- 
ing a value of the scattering length as an independent 
renormalization condition. This may appear too restric- 
tive and in fact it is possible to renormalize using energy 
dependent boundary conditions, a procedure essentially 
equivalent to imposing more renormalization conditions 
or counter-terms. Although there are subtleties on how 
short distances should be parameterized in such way that 
the cut-off may be removed [34|, |4l[ the procedure in co- 
ordinate space turns out to be rather simple. In the case 
of two conditions we would fix the scattering length, ao, 
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and the effective range ro independently of the poten- 
tial. The coordinate space procedure [34], [4j| consists of 
expanding the wave function in powers of the energy 

u p {r) = u {r) + p 2 u 2 (r) + ... , (44) 

where «o(r) and U2{r) satisfy the following equations, 



<(r) +MV(r)u (r) 



(45) 



uo{r) 



0, 

-> 1 - r/a , 

-u 2 \r) + MV{r)u 2 (r) = u (r) , (46) 
U2(r) — ► (r 3 - 3a r 2 + 3a r r) / (6a ) , 



The asymptotic conditions correspond to fix ao and ro 
as independent parameters (two counter-terms). The 
matching condition at the boundary r = r c becomes en- 
ergy dependent [34[ 



u ' P ( r c) _ u' (r c ) +p 2 u 2 {r c ) + 



u p (r c ) u (r c ) + p 2 U2(r c ) + 



(47) 



whence the corresponding phase shift may be deduced 
by integrating in Eq. (|4"5|) and Eq. ([4^]) and integrating 
out the finite energy equation. It is worth mentioning 
that the energy dependent matching condition, Eq. (|47|) . 
is quite unique since this is the only representation guar- 
anteeing the existence of the limit r c — > for singular 
potentials [34J- In any case, if ro is fixed from the start 
to their experimental values in the singlet and triplet 
channels, the Wigner correlation given by Eq. (f3"2")l and 
generating the universal curve shown in Fig. [T] would 
not be predicted and the symmetry between the 1 Sq and 
the 3 Si channels would be further hidden into the phase 
shifts. Note that the breaking of the symmetry with two 
counter-terms is a short distance one when the cut-off is 
eliminated, r c — > 0, since at any rate the potential is kept 
fixed and V s {r) — Vt(r) for any non- vanishing distance, 
r > r r > 0. Thus, if we write 



r = A + 



B 



C 



+ r 



;hort 



(48) 



with 



„short 



the effect of the second counter-term, we 



would obtain 



short 



short 



B 



1 

a, 



1 

a., 



1/a 2 terms have been neglected. 

„short 



(49) 



This 



where small 

yields rf hort — r s s hon ~ O.lfm. Thus, while introduc- 
ing no counter-term (trivial boundary condition) does 
not break the symmetry yielding identical phase shifts, 
5 s (k) = St(k), introducing more than one counter- term 
(energy dependent boundary condition) breaks the sym- 
metry at the ~ 10% level. As a consequence, we stick 
to the case of just one counter-term (energy independent 
boundary condition). 



B. Symmetry breaking due to the tensor force 

Of course, an interesting possibility which should be 
explored further is that of keeping the energy indepen- 
dence of the boundary condition and breaking the sym- 
metry by introducing a long distance component of the 
potential, such as e.g. the tensor force, which would in- 
clude the coupling of the 3 Si wave treated here to the 
3 Di channel. Actually, this would correspond to take 
into account, as proposed in Ref. [l6|, the leading and 
complete large-Ac NN potential. In other words, while 
Wigner symmetry implies a vanishing tensor force, lead- 
ing large- A c does not necessarily implies the tensor force 
to be small. To analyze this potential source of con- 
flict we consider the 3 Si effective range parameter which 
incorporates a D-wave contribution stemming from S-D 
tensor force mixing and is given by 



n = 2 



1 -u ,a(r) -w , a (ry 

at' 



dr . 



(50) 



where the zero energy S-wave function UQ a [r) — > 
u o, 3 Si( r ) (discussed above) and the D-wave function 
w o,a( r ) — > when the tensor force is switched off keep- 
ing a t fixed. The corresponding tensor potential would 
include tt and p exchange contributions characterized by 
the <7-7r n n and f P NN couplings and diverges as 1/r 3 at 
short distances. This situation resembles a previous OPE 
study 14611 and a detailed account will be presented else- 
where [47j. There, it will be shown how the extension 
of the superposition principle and renormalization to the 
coupled channel case yields in fact an identical analytical 
result as shown in Eq. (|3"2")l for the triplet (un-coupled) 
channel in the absence of tensor force. We will just 
quote here the numerical modification of the correlation 
relation coefficients for the triplet channel (the singlet 
1 So is not modified), Eq. (f3"2"| . Numerically, we get for 
fpNN = 17 and g^NN = 9-86 



r t = 2.6199 



5.7843 5.7608 



(51) 



which corresponds to a ~ 10% breaking due to the tensor 
force. As we see, the coefficients in Eq. (|33l) are not mod- 
ified much despite the singularity of the tensor force and 
its dominance at short distances. Actually, the depen- 
dence of the coefficients on the couplings responsible for 
the tensor force is moderate in a wide range. Therefore, 
while from the large N c viewpoint a large tensor force 
is not forbidden, we find the effect in the S-wave to be 
numerically small, as implied by Wigner symmetry. 

In this regard it should be noted that a virtue of the 
renormalization approach is that, since the scattering 
lengths are always fixed, such a long distance symmetry 
breaking term only influences the region where the poten- 
tial is resolved, and from this viewpoint the perturbation 
will be stable, i.e. the change will be small. Actually, 
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in Ref. [46| a suitable form of perturbation theory in the states, denoted as 5°j} and <5°°, holds, 
tensor force was suggested based on the known smallness 

of the mixing angle ei, which stays below 2 — 3°, in a wide iC^ + n h 

energy range and is indeed smaller than the 5p phase. It S L = — ■ — ^- - — - ^ = <5 iL = <5 L 

would be interesting to work out the consequences of such 



(21- + 1)3 

an approach when also p exchange is incorporated. „ n _ Y1j=l-i(?J + _ x oo _ roo /r / i\ 

^ L = (2L + 1)3 = hL = h ' (M) 



C. Symmetry breaking in non-central waves 

With the previous appealing interpretation of the 
Wigner symmetry as a long distance one for the S-waves, 
we analyze what are the consequences for the phase shifts 
corresponding to partial waves at angular momentum 
larger than zero, L > 0. Unlike the S-waves we expect 
the dependence on the short distance behavior to be sup- 
pressed due to the centrifugal barrier, and the symmetry 
should become more evident. Note also that while a dis- 
similarity between phase shifts connected by the symme- 
try does not necessarily imply long distance symmetry 
breaking, an identity between phase shifts is a clear hint 
of the symmetry. 

In the two-nucleon system the Wigner symmetry im- 
plies the following relations for spin-isospin components 
of the antisymmetric sextet, 6,4, and the symmetric decu- 
plet, 10s, respectively (see Appendix |A"|) thus we should 
have 

5°lj = SZ = S L , even-L (52) 
= Sfj = S L , odd - L (53) 

For P-waves, for instance, we have the spin singlet state 
1 P\ and the spin triplets 3 Po, 3 Pi and 3 p2 which accord- 
ing to the symmetry should be degenerate as they belong 
to the 10s supermultiplet. Inspection of the Nijmegen 
analysis [32| reveals that X P\ is very similar to 3 Pi at all 
energies, \Sip 1 — fopj ~ 1°, but very different from the 
3 Po and 3 p2 phases. For D-waves, associated to a 6^ su- 
permultiplet, we have a similarity between 1 £>2 and 3 Ds 
phases \Si D . 2 — 5a Da \ ~ 1° but, again, clear differences 
between the 3 Pi and 3 P/2 ones. Clearly, the symmetry is 
broken in higher partial waves. In what follows we want 
to determine whether our interpretation of a long dis- 
tance symmetry which worked so successfully for S-waves 
above (see Sect. IIV|) holds also for non-central phases. 

As it is well-known the spin-orbit interaction lifts the 
independence on the total angular momentum, via the 
operator L ■ S. Moreover, the tensor coupling operator, 
$12, mixes states with different orbital angular momen- 
tum. We proceed in first order perturbation theory, by 
using the Wigner symmetric distorted waves as the un- 
perturbed states. In appendix iDl wc show this procedure 
explicitly. To first order in spin-orbit and tensor force 
perturbation the following sum rule for the center of the 
5 = 1 multiplet, denoted by 8)® and 5}}, and the S = 



In terms of these mean phases, Wigner symmetry is for- 
mulated for non-central waves as 



Si Pl 


= l (&* + 


3c)3 Pl + 56a p 2 ) 


(55) 


Sl D 2 


= h 


+ 5^ + 7^3) 


(56) 


Si F3 


= k c 5 ** 


+ 7^3 +9<53 F J 


(57) 


^G 4 


= h {763G > 


+ 96 3Gi + ll&o 5 ) . 


(58) 



These sum rules are true as long as the short distance 
breaking can be considered small, and for this reason we 
have not written down the sum rule for S-waves. Fur- 
ther, they hold also when the tensor force is added. In 
Fig. [5] we show the l.h.s. and the r.h.s. of P-, D-, F- 
and G- waves. As we see the D-waves fulfill this relation 
rather accurately up to p ~ 250MeV and the G- waves up 
to p ~ 400MeV while the P- and F- waves fail completely. 
Actually, at threshold, Sl — > — a_Lp 2i+1 , and using the 
low energy parameters of the Nijmll and Reid93 poten- 
tials [12] determined in Ref. [48[ we get 

1 

m Pl = - [as Po + 3aa Pl + baspj , 
(~2.46fm 3 ) (0.08fm 3 ) , 

ai ° 2 = 15 ( 3a3£>1 + 5a3 £>2 + 7a3 D 3 ) , 
(-1.38fm 5 ) (-1.23ml 3 ), (59) 

where the numerical values are displayed below the sum 
rules. On the light of the previous discussions for the 
S-waves one reason for the discrepancy should be looked 
in a short distance breaking of the symmetry for the D- 
waves. Actually, the fact that D-waves violate the sum 
rule at p ~ 250MeV while the G-waves show no violation 
up to p ~ 400MeV agrees with our interpretation in the 
S-waves that the Wigner symmetry be a long distance 
one, since higher partial waves are less sensitive to short 
distance effects. The case of P-waves is different since 
the 1 P-potential and the 3 P-potentials are very different. 
This pattern of symmetry breaking agrees with the find- 
ings of Ref. [3] based on the large N c expansion where 
the central potential preserves the symmetry in L-even 
partial waves while it breaks the symmetry in the L-odd 
partial waves, since at leading order and neglecting the 
tensor force 

V(r) = V c {r) + arW s (r) + 0(1/N C ) , (60) 
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F-wave relation G-wave relation 




P cm [MeV] p cm [MeV] 



FIG. 6: Average values of the phase shifts [32I ] (in degrees) as a function of the CM momentum (in MeV) based on first order 
spin-orbit coupling. (Upper left panel) P-waves. (Upper right panel) D-waves. (Lower left panel) F-waves. (Lower right panel) 
G-waves. According to the Wigner symmetry 81 L — 5s L . Serber symmetry implies S3 L = for odd-L. One sees that L-even 
waves satisfy Wigner symmetry while L-odd waves satisfy Serber symmetry. 



so that for the lower L-channels we have 

Vi s = Vs s , = V c (r)-3W s (r)+0(l/N c ), 
Vi P = V c (r) + 9W s {r)+0(l/N c ), 
V. P = V c (r) + W s (r)+0(1/N c ), 

Vi D = V» D = V c {r)-3W s (r) + 0(l/N c ), 

(61) 

so as we see Vap ^ Vip, and thus it is obvious that 
S3 P =/= 5i P . One might check this further by proceeding 
as follows. In the case of odd waves such as the P-waves 
the proper comparison might be taking the 3 P-potential 
and renormalizing with the 3 P-mean scattering length, 
atap = 0.08fm 3 and compare to the 3 P-mean phase shift. 

We note that the initial claim of Ref. [l5| on the va- 
lidity of the Wigner symmetry based on the large N c ex- 
pansion was restricted to purely center potentials, which 
do not faithfully distinguish the two irreducible repre- 
sentations, IO5 and 6a, of the 5f (4) group for the NN 



system. Later on, the issue was qualified by a more com- 
plete study carried out in Ref. [Ty| which in fact could 
not justify the Wigner symmetry in odd-L partial waves, 
even when the tensor force was neglected. Although this 
appeared as a puzzling result, it is amazing to note that 
our calculations clearly show that the pattern of SU (4)- 
symmetry breaking supports a weak violation in even-L 
partial waves and a strong violation in the odd-L partial 
waves, exactly as the large N c expansion suggests. 



D. Serber symmetry 

On the other hand, from the odd-waves we see from 
Fig.[5]that the mean triplet phase is close to null, thus one 
might attribute this feature to an accidental symmetry 
where the odd-waves potentials are likewise negligible. In 
the large N c limit this means Vc + Ws 3> Vc + 9Ws, a 
fact which is well verified. For instance at short distances 
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the Yukawa OBE potentials have Coulomb like behavior 
V — > C/(47rr) with the dimensionless combinations 



Cvc+Ws 



2 2 
~9aNN + 5ojATAT + 



6M2 



Vc+9Ws 



2 2 
~9aNN + 9uNN + 



2M N 



(62) 

where the small OPE contribution has been dropped. 
Numerically we get Cv c +w s ~ 10 and Cv c +9W S ~ 300 
for reasonable choice of couplings. Although this approx- 
imate vanishing of triplet odd-wave potentials is not a 
consequence of large iV c it is nevertheless reminiscent of 
the old and well-known Serber force, 

^Serber(r) - ^ (1 + P M ) \ (1 - P a ) V s (f) 

+ \{t-PM)\{l + P„)V t {r), (63) 

with Pm the Majorana coordinate exchange operator. 
Due to the Pauli principle PmPoPt = —1 with P T = 
(1 + t)/2 and P a = (1 + er)/2 the isospin and spin ex- 
change yields vanishing potentials for spin-triplet and 
isospin-triplet channels, and generating a scattering am- 
plitude which is even in the CM scattering angle, a prop- 
erty which is approximately well fulfilled experimentally 
for pp-scattering. We call this property Serber symmetry 
for definiteness. After introducing spin-orbit coupling we 
would get the sum rules to first order 



3p 



9 



(<53 Po + 3* 



3 Pi 



5(5; 



■■P-2 



= 0. 



Si 



<F 
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i F2 + 7d3 F3 + 9d3 Fi 



0. 



(64) 
(65) 



which is well fulfilled by the phase shifts [32| as shown in 
Fig. [6] where 63 P <C dip and 5a F <C Si F . In the large N c 
limit we may comply both with Wigner symmetry in L- 
even waves and Serber symmetry in L-odd waves when 
Ws(r) — —Vc(r), whence generally V(r) — Vc(r)(l — 
err). Even if we neglect the small OPE effects, this will 
clearly not be exactly fulfilled unless one would require 
vn p = = m CT . Although there are schemes where such 
an identity between scalar and vector meson masses are 
explicitly verified [H, HH, at present, it is unclear 
whether the Serber symmetry which we observe in the 
NN system for spin-triplet and odd-L phase-shifts could 
be formulated as a symmetry from the underlying QCD 
Lagrangean. 

Our findings suggest that a pure large N c in the ab- 
sence of tensor force not only is compatible with the 
standard Wigner symmetry in the case of the dominant 
S- waves and higher L-even channels, but it might also 
be a competitive alternative for the L-odd waves where 
the usual Wigner symmetry is broken and Serber sym- 
metry holds instead. Of course it would be interesting 
to pursue the more complete situation including the ten- 
sor force from the start, a case which will be presented 
elsewhere 47 1. 



E. NN Level density in the continuum 

Our results have some impact for hot nuclear matter 
at low densities. In the continuum, we may think of 
putting the two nucleon system in a box and evaluate 
the corresponding level density when the infinite volume 
limit is taken. This is a standard problem in statistical 
mechanics which appears, e.g. in the calculation of the 
second virial coefficient contribution to the equation of 
state of a dilute quantum gas [13] (see Refs. [H, [54| for 
recent applications to hot nuclear matter). The result is 
expressed as 



p{E) = -Li* 
2m 



S(E) 



dE 



1 dA NN (E) 
7T dE 



(66) 



where S(E) is the S matrix in all coupled channels and 
the total phase A is defined by 

A NN (E) = ( 2 J + !)( 2T + Wlj(Z) ■ (67) 

S,T,J 

In the case of coupled channels one should consider the 
corresponding cigenphases 12 . Defining the mean phase 
as 



Ejt S L -s^J+^l T j(E) 
(25+l)(2L + l) 



5l T {E) = 



(68) 



corresponding to the phase-shift analog of the center of 
gravity of the supermultiplet (see also Eq. ([54]) we get 

A NN (E) =Y,( 2S + !)( 2L + X )( 2T + !)^ T (^) ■ (69) 

S,T,J 

Thus, using the above relations, Eq. (|58p for L-even waves 
and Eq. (|6"5]) for L-odd waves, featuring Wigner and Ser- 
ber symmetries respectively we would get that mixed 
triplet channel contributions may be either eliminated 
in terms of singlet ones for even-L or do not contribute 
for odd-L, 

A NN (E) = 3 {Si So + Ss Sl ) + 3Si Pl + 30Si D2 +... (70) 

For the neutron case we have 

A nn (E) = 5i Sq + 5Si D2 + 98i G4 + ... (71) 

i.e., odd-L waves do not contribute. The lack of a P-wave 
contribution scaling as ~ —ctpp 3 is compatible with the 
minimum observed in Ref. [53| for A„„ in the subthresh- 
old region E'lab < 50MeV. 



In the special case of NN scattering one can also use the nu- 
clear bar phase shifts due to the identity 



S3 



(J 



neglecting 



)j 



(■74 



j\ . The concern spelled out in Ref. [53(| that 
the mixing was an approximation is unjustified. 
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VI. CONCLUSIONS 

At low energies NN interactions are dominated by two 
S-waves in different channels where spin-isospin (S, T) 
are interchanged, (1,0) <-> (0,1). Wigner SU(A) symme- 
try implies that the potentials in the 1 So and 3 Si chan- 
nels coincide and the tensor force vanishes, while the cor- 
responding phase shifts from Partial Wave Analyses are 
actually very different at all energies and show no evident 
trace of the identity of the potential, besides the qualita- 
tive fact that a weakly bound deuteron 3 Si state and an 
almost bound virtual So take place. Given the fact that 
the nuclear force at short distances is fairly unknown, 
the validity of the symmetry to all distances would be at 
least questionable and could hardly be tested quantita- 
tively. On the other hand, our lack of knowledge of the 
short distance physics should not be crucial at low ener- 
gies, where the phase shifts are indeed quite dissimilar. 
Therefore, we propose to regard St/(4) as a long dis- 
tance symmetry which might be strongly broken at short 
distances and weakly broken at large distances. Using 
renormalization ideas where the desirable short distance 
insensitivity is manifestly fulfilled we have shown how the 
standard Wigner correlation between potentials indeed 
predicts one phase shift from the other in a non-trivial 
and successful way. Remarkably, using a large N c moti- 
vated One Boson Exchange potential we have proven that 
if one channel is described successfully the other chan- 
nel is unavoidably well reproduced within uncertainties 
which might be compatible with the disregard of the ten- 
sor force and the 1/N% corrections to the potential. This 
long distance correlation holds also for the virtual singlet 
state and the deuteron bound state. Actually, the effects 
of symmetry breaking at long and short distances have 
been analyzed and the extension to higher partial waves 
has also been discussed, where a relation for phase shifts 
has been deduced. Our calculations provide a justifica- 
tion on the use of Wigner symmetry in light nuclei solely 
on the basis of the NN-interaction and suggest that a 
specific interpretation of the Wigner symmetry as a long 
distance one in conjunction with renormalization theory 
extends beyond the scaling region to a much wider range 
than assumed hitherto. It would be interesting to see how 
these ideas could be further exploited beyond the simple 
two nucleon system. However, a key question has always 
been what is the origin of the accidental Wigner symme- 
try from the underlying fundamental QCD Lagrangean 
and, moreover, under what conditions this is expected to 
be a useful symmetry. We find the large N c expansion 
in the absence of tensor force besides being compatible 
with the standard Wigner symmetry in the case of the 
low energy dominant S-waves and subdominant higher 
L-even partial waves it may also become a competitive 
alternative for the other L-odd partial waves where the 
usual Wigner symmetry is manifestly broken. These con- 
clusions are remarkable, for they suggest that a unfore- 
seen handle on the nature, applicability and interpreta- 
tion of a widely used approximate nuclear symmetry may 



be based on a QCD distinct pattern such as the large N c 
limit. Obviously, it would be very interesting to pursue 
further the study of the complete large N c potential with 
inclusion of the tensor force to verify this issue in more 
detail [47] . In our view this would definitely provide use- 
ful insights into QCD inspired approximation schemes in 
nuclear physics. 
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APPENDIX A: WIGNER SYMMETRY FOR NN 



Wigner ST/ (4) spin-isospin symmetry consists of the 
following 15-generators 0,0, EJ 



Z A 
L A 

Gia \ ^ a 



(Al) 
(A2) 
(A3) 



where t\ and a\ are isospin and spin Pauli matrices for 
nucleon A respectively, and T a is the total isospin, S l the 
total spin and G la the Gamow- Teller transition operator. 
The quadratic Casimir operator reads 



C 



SU{i) 



T a T a + + G ia G l: 



(A4) 



and a complete set of commuting operators can be taken 
to be CsuU)t T3 and S 2 ,G 2 3. The fundamental repre- 
sentation has Csu(4) — 4 and corresponds to a single 
nucleon state with a quartet of states p ] , p n ] , n {, 
with total spin S —1/2 and isospin T =1/2 represented 
4 = (S,T) = (1/2,1/2). For two nucleon states with 
good spin S and good isospin T Pauli principle requires 
(— ) S+T+L — — 1 with L the angular momentum, thus 



ST 



1 



SU(4) = g ( CT + T + CTT ) + — • 
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y 



(A5) 



where r = t\ ■ t% — 2T(T + 1) — 3 and a = a± ■ &i — 
2S(S -I- 1) — 3 and the corresponding wave function is of 
the form 



,LS 



(r) 



Y, 



LM 



,{i)x SMs x™ T 



(A6) 
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One has two supermultiplets, which Casimir values are 



/-<oo _ /'-ill _ o 

°S£/(4) — °SC/(4) — y ) 



/-lOl 
°S£/(4) 



/-ilO 
°S£/(4) 



5, 



(A7) 
(A8) 



corresponding to an antisymmetric sextet 6a — (0, 1) 
(1,0) when L = even and a symmetric decuplet 10g = 
(0,0) (1, 1) when L =odd. The radial wave functions 
fulfill ?i" 1 (r) = v}®{r) and it°°(r) = uj}(r) respectively. 
This means that we have the following supermultiplets 

CS , 3 Si) , CP!, 3 P ,i j2 ) , CD 2 , 3 £> 1>2 , 3 ) . . . (A9) 

When applied to the NN potential, the requirement of 
Wigner symmetry for all states, implies 

V T = Wt = V ls = W ls = 0, 
W s = Vs = W c , (A10) 



so that the potential may be written as 



V = V c + (2Cf£ (4) - 15)W S 



(All) 



Note that the particular choice Ws = corresponds 
to a spin-isospin independent potential, but in this case 
no distinction between the 6a and IO5 supermultiplets 
arises. As it is well-known the spin-orbit interaction lifts 
the total angular momentum independence. The Wigner 
symmetry does not distinguish between different total an- 
gular momentum values, so admitting that the potentials 
are different we may define a common potential 



V LST {r) 



sr^L+S 

2-ij 



L-S 



(2J+l)V JST (r) 



(25+l)(2i + l) 



where similarly to the perturbation theory for energy lev- 
els where the center of a multiplet of states is predicted, 
the appropriate statistical weights related to the angu- 
lar momentum have been used. The previous expression 
makes sense if the symmetry is broken linearly by spin- 
orbit coupling. In terms of these mean potentials the 
symmetry would be 



or equivalently 



Vr L (r) = V3 L (r), 



3(2L+ 1) 



(A13) 



(A14) 



As mentioned in the paper, if the symmetry is taken lit- 
erally at all distances we should have 5il — 8^ L . 



APPENDIX B: LONG DISTANCE 
PERTURBATION THEORY 

We illustrate here a situation where the potential may 
be treated in long distance perturbation theory and 



renormalized (for a somewhat similar approach for finite 
cut-offs see e.g. Ref. [55]). Unlike the standard per- 
turbative approach, which usually does not hold in the 
presence of bound states, this expansion can deal with 
weakly bound states, provided this is the only one. This 
is in fact the case for the OPE potential for the parame- 
ters we use, applied to the deuteron state, for which we 
show the procedure here to first order. To analyze this 
situation we vary the potential V — > V + AV 



Au k (r)" 



MAV(r)u k (r) 
MV(r)Au k (r) 



k 2 Au k (r), (Bl) 



we use the previous wave functions u k (r) as the zeroth 
order approximation, corresponding to take V(r) = 
and solve for the first order correction Au k (r) the equa- 
tion which asymptotic wave function corresponds to take 
the phase shift 8 + AS. Multiplying Eq. © by Au k (r) 
and Eq. (|B1[) by u k (r), subtracting both equations and 
integrating from r c to 00 we get 



i' k Au k + u k Au k ] 



dr AU{r)u k {r) 2 . (B2) 



The lower limit term may be related to the variation of 
the boundary condition, whereas the upper limit term 
is related to the change in the phase shift, AS. In or- 
der to eliminate the cut-off we subtract the zero energy 
limit, k — > 0, and using the energy independence of the 
boundary condition we get some cancellation since 



u k (r c ) u [r c ) 



(B3) 



(A12) 

Finally, the result may be re-written as follows 



A (kcotS) = -A I — 



(B4) 



AU(r) \u k {rf - u (r) 2 ] dr . 



If we fix the scattering length independently on the po- 
tential we have Aao = 0, thus eliminating the first term 
of the r.h.s. and after taking the limit r c — > the result 
for the total (and renormalized) phase shift to first order 
in the potential reads 



1 r 

kcotSo(k) = h / 

ao Jo 



cos(fcr) 



drMV(r) 
sin(fcr) 



(B5) 



otok 



r 

on 



The renormalized effective range is entirely predicted 
from the potential at all distances 



ro 



drr 2 MV{r) ^1 - +... (B6) 



Note the extra power suppression at the origin when ao 
is fixed independently on the potential, indicating short 
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distances become less important. The bound state can be 
obtained in a similar manner by replacing u/-(r) — > u 7 (r), 
assuming that the binding energy is independent on the 
potential, A7 = 0, and using orthogonality Eq. (|B3[) to 
the zero energy state 

1 f°° 

— = 7 + / MV{r) \uJrf - u (r) 2 ] dr . (B7) 
«o Jo 

This equation is implicit in both ao and 7, but we can 
make it perturbative explicitly, using that to first order 
a o ~ I/7 in the zero energy wave function u (r) ~ 1 — jr, 
yielding 

1 r 00 

— = 7 + / MV(r) \e' 2 ^' - (1 - 7 r) 2 l dr (B8) 
"o Jo 



APPENDIX C: SCALE INVARIANCE AND 
RENORMALIZATION 

We have suggested that Wigner symmetry be a long 
distance one. From a renormalization group (RG) view- 
point this has a simple interpretation (for a discussion in 
coordinate space see e.g. Ref. [H], It means finding 
a solution to the RG equations which break the symme- 
try of the equations. A very simple case which illustrates 
this issue is provided by the problem 



k 2 u(r) . 



(CI) 



At zero energy, k = 0, the solution is invariant under the 
scaling transformation r — » Ar. This property holds also 
at short distances, where the energy term on the r.h.s. 
can be neglected. If we use the RG equation, Eq. ([35]) . 
for this particular case at short distances 



Rdo{R) = co(fl)(l - co(-R)) + g . 



(C2) 



The scale symmetry becomes now evident; if co(R) is a 
solution then cq(XR) is also a solution for any value of 
A 7^ 0. The solution must necessarily specify the value at 
a given scale co(Ro), hence breaking explicitly the dilata- 
tion symmetry. This symmetry breaking is unavoidable. 
In Refs. H] it is shown how, for g < —1/4 the break- 
ing is lowered to the discrete subgroup of dilatations, and 
the connection to the Russian Doll renormalization. In 
the case of the Wigner symmetry for the 1 S'o and 3 So 
potentials discussed in the paper, the breaking is not un- 
avoidable, and there exists in fact a very special choice 
where the symmetry can be preserved by taking identical 
boundary conditions at a given scale. Besides this par- 
ticular solution, the identity between solutions Co tS (R) 
and co,t(i?) will generally be violated, although the rela- 
tion from one scale to a different one co, s (-Ro) — > Cq, s (R) 
and co,t(i?o) ~* co,t(iJ) is governed by the same relation, 
Eq. (Jp). 

It is worth noting the resemblance of the previous 
quantum-mechanical discussion with similar and well- 
known field theoretical concepts. The unavoidable 



breaking of the dilatation symmetry corresponds to an 
anomaly of the dilatation current. The optional choice 
of boundary conditions corresponds to the case of finite 
but ambiguous theories (see e.g. Ref. |4pj). 



APPENDIX D: SPLITTING FORMULA FOR 
PHASE-SHIFTS 

We want to derive the splitting formula for phase shifts, 
Eq. (|54[) by using distorted waves perturbation theory. 
The coupled channel Schrodinger equation for the rela- 
tive motion reads 



u"(r) 



U(r) + -v 



u(r) = k 2 u(r) . 



(Dl) 



where Uf J L ,(r) = 2ii np Vf J L ,(r) is the coupled channel 
matrix potential which for the total angular momentum 
J > can be written as, 



U 0J (r) = uy } 



O.J 



(D2) 



U 1J (r) = 



fTTlJ 
U J-l,J 



x(r) 
UYj{r) 




TTlJ 

U J-1,J^ 



TTlJ 
U J+1,J- 



l( r ) N 

i(r) y 



In Eq. dDl} L 2 = diag(Li(£i + 1),...,L N (L N + 1)) is 
the angular momentum, u(r) is the reduced matrix wave 
function and k the CM. momentum. In the case at hand 
N = 1 for the spin singlet channel with L = J and N = 3 
for the spin triplet channel with L\ = J — 1, Li = J 
and L3 = J + 1 . For ease of notation we will keep the 
compact matrix notation of Eq. (|D1|) . At long distances, 
we assume the asymptotic normalization condition 



u(r) -> hH(r)-h(+)( r )S, 



(D3) 



with S the standard coupled channel unitary S-matrix. 
For the spin singlet state, S — 0, one has L = J and 
hence the state is un-coupled 



whereas for the spin triplet state S 
un-coupled L = J state 



-.1.7 



= e 2lS1 '' 



(D4) 

1, one has the 
(D5) 



and the two channel coupled states L, L' = j ± 1 states 
which written in terms of the eigenphases are 



S 1J = 



/ cos e j — sin ej 
I sine j cosej 



cos e j sin ej 
— sin e j cos e j 








(D6) 



The corresponding out-going and in-going free spherical 
waves are given by 



,(±: 



(r) = diag(/if i (fcr), 



,(fcr)), (D7) 
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with h^(x) the reduced Hankel functions of order I, 



hi 



(,:) 



L + l/2 



(x) ( ht 



and satisfy the free 



Schrodinger's equation for a free particle. 

In order to determine the infinitesimal change of the 
S matrix, S — > S + AS, under a general deforma- 
tion of the potential U(r) — > U(r) + AU(r) we use 
Schrodinger's equation (|Dlj) and the standard Lagrange's 
identity adapted to this particular case, we get 

[u(r) T Au'(r) - u'(r)t Au(r)] ' = u(r) t AU(r)u(r) . 

(D8) 

The unitarity of the S-matrix, S T S = 1, yields the con- 
dition AS^S + AS = 0. We assume a mixed boundary 
condition at short distances, r = r c , for the unperturbed 
coupled channel potential, U(r), 



u'(r c )+Lu(r c ) =0, 



(D9) 



with L a self-adjoint matrix. After integration from the 
cut-off radius r c to infinity and using the asymptotic form 
of the matrix wave function, Eq. (|D3|) . as well as the 
condition at the origin, Eq. (|D9j) yields 



/>oo 

2ifcS f AS= / dru(r) + AU(r)u(r). (D10) 

If we take the Wigner symmetric states as the unper- 
turbed problem, then S , U(r) and u(r) become a diag- 
onal matrices, so that 

ASS ' I = -^L ^f>)tAU(r)uf>), (Dll) 



so that the perturbed eigenphases become 



rST _ cST 



AS 



ST 
JL 



(D12) 



Note that to this order the mixing phases vanish, Aej = 
0. Identifying further AU with the spin-orbit and the 
tensor potential, in the LS-coupling the result may be 
written as 



8E 



5l T 
A S L T 



S,1^L 

[J(J + l)-L(L + l)-S(S + 



1)] , 
(D13) 



where (Sj^Jj-i.j-x = -2( J - l)/(2 J + 1), (Sf 2 )j,j = 2, 
(5/2)7+1,7+1 = -2(J+ 2)/(2J+ 1). Defining the super- 



multiplet coefficients Al 



Al 1 



Af 



A / and B L = A° T ° 



rlO 
°L,7 



A L [J(J + 1) - L(L + 1) - 2] 



Cl{S 12 )ll , 
5 L ° J + B L [J(J- 

Dl(S( 2 )ll , 



(D14) 



1) - L(L + 1) - 2] 



(D15) 



we readily get the sum rule for phase-shifts, Eq. (|54[) . 
The above equations would yield a Lande-like interval 
rule between spin-triplet energy levels for the spin-orbit 
or the tensor potentials separately. For instance, 



5i Pl 

63p 2 



Sp, 

8 P - 4£>i - ABx , 
S P + 2D 1 - 2Bx , 



S P - -Di 
5 



2Bi 



(D16) 



A further remark is in order, since the spin-orbit or tensor 
potentials may be singular at the origin. In such a case 
of singular perturbations one computes the sum rule first 
and then removes the cut-off, r c — ► 0. 
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